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ABSTRACT 

We present extensive spectral and photometric observations of the recycled pulsar/white-dwarf 
binary containing PSR J0437— 4715, which we analyzed together with archival X-ray and gamma- 
ray data, to obtain the complete mid-infrared to gamma-ray spectrum. We first fit each part of 
the spectrum separately, and then the whole multi-wavelength spectrum. We find that the optical- 
infrared part of the spectrum is well fit by a cool white dwarf atmosphere model with pure hydrogen 
composition. The model atmosphere (Teff = 3950 ± 150K, \ogg = 6.98 ± 0.15, i?wD = (1.9 ± 0.2) x 
10^ cm) fits our spectral data remarkably well for the known mass and distance (M = O.25±O.O2M0, 
d = 156.3 ± 1.3 pc), yielding the white dwarf age (twd = 6.0 ± 0.5 Gyr). In the UV, we find a 
spectral shape consistent with thermal emission from the bulk of the neutron star surface, with surface 
temperature between 1.25 x 10^ and 3.5 x 10^ K. The temprature of the thermal spectrum suggests 
that some heating mechanism operates throughout the life of the neutron star. The temperature 
distribution on the neutron star surface is non-uniform. In the X-rays, we confirm the presence 
of a high-energy tail which is consistent with a continuation of the cut-off power-law component 
(r = 1.56 ± 0.01, i?cut = 1-1 ± 0.2 GeV) that is seen in gamma-rays and perhaps even extends to the 
near-UV. 

Subject headings: pulsars: individual (PSR J0437— 4715); white dwarfs 



1. INTRODUCTION 

Millisecond pulsars (MSPs) are neutron stars that have 
been spun up through the accretion of matter from the 
binary companion at some time in the past, as part of a 
binary system. This process, known as recycling, enables 
the pulsar to emit detectable emission long after it would 
have died as a classical pulsar. 

MSPs have mainly been observed in the radio, where 
the emission is dominated by coherent processes, and 
in X-rays and 7-rays. The X-ray radiation is thought 

lmartin.dtu'ant®astro.ufl.edul 



to be comprised of two components: non-thermal emis- 
sion from the magnetosphere, which is characterized by 
hard spectra (photon index F = 1 — 2) and sharp pul- 
sations with high pulsed fraction; and thermal, surface 
emission, which has softer spectra and smoother pul- 
sations. The former mechanism dominates in MSPs 
with high spin-down luminositjQ E > lO'^^ergs"^ (e.g., 
PSR B1821-24 and B1937+21), whereas a second group 
of MSPs with E = 10^3 - lO^^ergs"! (e.g., PSR 

^ Rate of loss of rotational kinetic energy E = IQfl, for a stan- 
dard neutron star moment of inertia I = 10^^ gcm^. 
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J0030+0451) show mainly thermal emission, with a pos- 
sible power-law (PL) high-energy component. For ther- 
mal spectra, X-ray observations are most sensitive to 
temperatures T ^fewxlO^K, which are considered to 
be typical temperatures for MSP polar caps, heated by 
relativistic particles. 

PSR J0437-4715 is an MSP (P = 5.8 ms) in a bi- 
nary syste m with a w h ite dw arf (WD), with orbital pe- 
riod 5.5 d. iBell et all ()1995f ) measured the proper mo- 
tion of PSR J0437-4715, 135±4mas/yr and detected 
an Hq, bow-shoc k nebula in t he approximate direction 
of motion; later iDeller et al.i (|2008) greatly improved 
the accuracy of the proper motion measurements using 
VLBI {fic, = 121.679 ± 0.052 mas/yr, ^ls = -71.820 ± 
0.086 mas/yr), and provided the most accurate distance 
to PSR J0437-4715, d = 156.3±1.3pc (making it one of 
the closest pulsars known). The masses of the two bi- 
nary components are Mpsr = 1.76 ± .20 Mfr^; Mwd = 
0.254 ± 0.018 Mq (jVerbiest et al.l [20081) . and the orbital 
inclination i = 137.6 ± 0.2°. PSR J0437-4715 is the 
brightest of the MSPs in the UV and X-rays. Detailed 
spectral studies have suggested that the bulk of the X-ray 
emission comes from hot polar regions, where the tem- 
perature is no n-uniform, decreasing with distance awa y 
from the pole (jZavlin fc Pavlovlll998HZavlin et al.ll200l . 

For the WD companion. iDanziger et al" I (|1993n~ give 
black-body fit of T '^4000 K to ground-based broad-band 
photometry. With such a cool WD, its emission no longer 
dominates in the UV if the NS is hot enough. This offers 
a unique opportunity to observe thermal emission from 
t he bulk of the neutron s tar surface in the UV/optical. 

iKargaltsev et all ()2004[ ) (hereafter K04) observed PSR 
J0437-4715 for the first time in the UV with the Space 
Telescope Imaging Spectrometer (STIS) aboard HST. 
K04 found a spectrum consistent with the Rayleigh- Jeans 
tail of thermal emission although the measured spectral 
shape had a large uncertainty. Assuming the maximum 
area that this could be emitted from is the surface of 
the neutron star, they derived a surface temperature of 
10^ K. 

Here we describe an extensive observational campaign 
ofPSR J0437-4715 in various spectral bands. IN Sec- 
tion 2 we descibe the data and its analysis, one spectral 
window at a time. The results and spectral fits from 
each dataset are presented in Section 3, with a complete 
multi-wavelength analysis in Section 4. The results are 
discussed in Section 5, followed by a brief summary and 
outlook in Section 6. 

2. OBSERVATIONS AND ANALYSIS 

We utilize a suite of observations of PSR J0437— 4715 
from ground- and space-based observatories. Here we de- 
scribe in detail each part of our campaign, and the qual- 
ity of the resulting data. We obtain spectral data in the 
mid-IR to 7-rays, with a combination of spectroscopic 
and imaging photometry observations (see Table 1). We 
made use of ACS on HST, FORSl on VLT/UTl of 
ESO/Paranal, and PANIC on Magellan I. Furthermore, 
we have re-analyzed archival Chandra, XMM-Newton, 
Spitzer and Fermi data. The only measurements al- 
ready in the literature which we have included (without 
reanalysis) are the gr ound-based photometric points of 
IDanziger et all (pM) . and the STIS FUV fluxes of K04. 



TABLE 1 

Log of observations of PSR J0437— 4715 



Waveband 


Observatory/ 
Instrument 


Date 


Type 


Section 


mid-IR 


Spitzer/MIPS 


Sep 2005 


Phot 


2.1.2 


mid-IR 


Spitzer/IRAC 


Nov 2004 


Phot 


2.1.1 


near-IR 


Magellan/PANIC 


Oct 2005 


Phot 


2.2 


optical 


VLT/FORSl 


Jul 1999 


Spec 


2.2 


optical 


HST/WFPC2 


May 1996 


Phot 


2.4.1 


op/UV 


HST/ACS 


Apr 2006 


Phot 


2.4.2 


NUV 


HST/ACS 


Apr 2006 


Spec 


2.4.3 


FUV 


HST/ACS 


May 2006 


Spec 


2.4.4 


X-ray 


XMM/E-MOS 


Oct 2002 


ImSpec'' 


2.5.1 


X-ray 


Chandra/ ACIS 


various 


ImSpec^ 


2.5.2 


Gamma 


FERMI/LAT 


2008-2010 


ImSpec^ 


2.6 



^ imaging on a detector with intrinsic energy resolution. 

2.1. Spitzer mid-IR photometry 

There are archival IRAC and MIPS Spitzer observa- 
tions of PSR J0437— 4715. The observations in the four 
Spitzer IRAC cha nnels, 3.6 /um, 4.5 /xm, 5.8^m, 8.0 ^m, 
(jFazio et al.|[200l were carried out in November 2004, 
the MIPS 24 /xm observations (jRieke et all I2004D were 
done in September 2005. The pointing accuracy of 
Spitzer is reported to be better than 0'.'5, according to 
the IRAC instrument handbool0. Astrometric accuracy 
of better than 0'.'2 can be reached, especially if cross- 
correlating with 2MASS point sources is successful. This 
is the case for the observations of PSR J0437— 4715, 
where around 20 2MASS point sources with IRAC coun- 
terparts are in the field. We apply the Spitzer Software, 
MOPEX (version 18.3.30 and its point source extrac- 
tion package, APEX, for the data reduction, which we 
describe in detail. 

2.1.1. IRAC observations 

The IRAC observations are strongly affected by ex- 
tended emission, which turns out to be an extended arti- 
fact in the detector, probably latent from previous obser- 
vations. We select only those Corrected Basic Calibrated 
Data files (CBCD; they already include some pipeline 
processing) from the dithered observations where the tar- 
get position is well separated from the artifact. This re- 
duces the depth of the observations; in the case of the 
fourth channel (8 /im), for example, only half of the avail- 
able coverage can be used. For this channel we removed 
one more CBCD file because of a cosmic ray hit right at 
the target position. The resulting mosaics at the target 
position are shown in Figure [1] for the first channel the 
position at the epoch of the IRAC observation of PSR 
J0437-4715 is indicated. 

Spitzer IRAC photometric accuracy can reach 2% for 
bright sources, provided the proper aperture correc- 
tion, the correction of the photometry array location 
effect, and the col or correction are taken into account 
(jReach et al.l I2005() . As recommended by the Spitzer 
Science Center, we use aperture photometry instead of 
profile fitting photometry for the measurements, using a 
source aperture radius of 2 native IRAC pixels (1'.'2 pixel 
size) and a sky annulus from 2 to 6 native IRAC pix- 
els to avoid bright sources in the neighborhood. (This 

http:/ /ssc. spitzer. caltech.edu/irac/iracinstrumenthandbook 
^ http:/ /ssc. spitzer. caltech.edu/postbcd/mopex. html 




Fig. 1. — IRAC images of the field of J 437- 4715 in four channels. The images are 36 X 30 in size. The VLBI radio position of 
PSR J0437— 4715 obtained by IDeller et al.l II2008I ). and corrected for the proper motion to the epoch of the Spitzer IRAC observations, 
is marked by a white cross in IRAC channel 1. The solid circle indicates the aperture (radius of 2 native IRAC pixels) used for the flux 
density measurements. The dashed circle corresponds to the aperture measurement done for the neighboring source Nl in channel 2. The 
same position and aperture was also used for channel 1 (see text), while in channels 3 and 4 Nl is not seen. 



results in a slight over-estimation of the sky due to flux 
from the source in the sky aperture. For the bright back- 
ground and close neighboring sources, the effect is small 
compared to the photometric uncertainties.) We use the 
'mode' algorithm for the background estimation in the 
annulus to limit the influence of a faint neighbor cover- 
ing part of the annulus. The APEX source extraction 
clearly detects a source at the target position in IRAC 
channels 1 and 2, the source is marginally detected in 
channel 3, and not detected in channel 4. The following 
flux values for channels 1 to 3 are aperture corrected 
as well as corrected for the position-dependent sensi- 
tivity, following the IRAC data handbook recommenda- 
tions: 16.0 ± 0.4 ^Jy; 11.4 ± 0.7 ^fJy; and 8.2 ±2.5 ^fJy for 
channels 1, 2 and 3, respectively, where uncertainties are 
at the Icr-level. The color correction for a 4000-5000 K 
blackbody is 1% or less for each channel, so we neglect 



this. 

The choice of the aperture and annulus size has been 
optimized to include as few neighboring stars as pos- 
sible. At larger separations from the target additional 
constraints for background estimates are the remaining 
artifacts and the varying IR background (we found back- 
ground differences of up to a factor of 3). For IRAC 
channel 1 we used apertures and annuli of the same size 
at three different positions in the immediate vicinity of 
the target to check for conspicuous deviations in the es- 
timated background or noise estimates. The obtained 
background and noise values are comparable to those of 
the target measurements. Therefore, we deduce that our 
aperture and annulus choice is acceptable with no ex- 
pected significant systematic effects regarding the back- 
ground or noise estimates. For the longer wavelengths, 
the variation in the IR background is expected to be even 
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higher. 

In IRAC channels 1 and 2 there is a neighboring source 
close to the target. This source is not seen in channels 
3 or 4. The applied apertures for the target as well as 
for this other source are shown in Figure [TJ in case of 
channels 3 or 4 no aperture measurements of the neigh- 
bor source were required, and the dashed circles are only 
plotted to indicate its position. The APEX source ex- 
traction detects the source only in channel 2; we use this 
position to do aperture photometry in channel 1 as well. 
The flux density of the neighbor source in channel 2 is 
8.5 ± 1.0 /zJy and for channel 1 is 6.2 ± 1.2 /iJy (Icr con- 
fidence) . 

Roughly 1 /3 of the neighbor aperture area is also cov- 
ered by the target aperture. If the apparent neighbor 
were a real source, the flux of our target would be the 
measured value minus the neighbor's contribution; if it 
were instead some manner of leakage, the target flux 
would be the sum of both fluxes. As we don't know 
which of the two options takes place, we leave the flux 
unchanged, and increase the uncertainty by 1/3 of the 
apparent neighbor flux. Thus, the target flux density in 
channel 2 is 11.4±3.3/iJy. Following the same argument 
as for channel 2, the target flux density in channel 1 is 
16.0±2.7/LtJy. 

We note that the neighboring source is not seen in any 
of our other observations (X-ray to NIR). The 3f7 Kg- 
band limit is 20.7 mag (see below), indicating an emis- 
sion maximum at IRAC wavelengths for this object. We 
suspect, therefore, that it is not a physical sourcfl 

The target is not detected in channel 4. We use the 
source position from IRAC channel 1 for aperture pho- 
tometry in order to obtain an upper limit for channel 4. 
We derive a 3cr limit of 22.6 /iJy on the source flux. 

2.1.2. MIPS 24 fim observations 

The photometric properties of MIPS 24/im observa- 
tions were reported to be excellent with achievable RMS 
errors of 2% and an accura cy of the absolute calibration 
being better than 5% tRieke et all [2004: MIPS instru- 
ment handbool^B)- The pipeline processed PBCD (post- 
BCD) 24 fj.m. mosaic of PSR J0437-4715 does not show 
any obvious artifact at the position of the target. As 
a check, we self-calibrated the data by dividing each 
EBCD (enhanced BCD) image with a normalized flat 
field created by taking the median of all the EBCDs. 
Self-calibration in this way removes large-scale sky vari- 
ations and some artifacts. In our case, self-calibration did 
not improve the data, and we use the unchanged PBCD 
mosaic from the Spitzer pipeline (version: S18.12.0) for 
our analysis. At the position of the target there seems 
to be slightly increased count rate; however, it is close to 
the noise level. We use MOPEX multi-frame point re- 
sponse function (PRF)-fitting photometry and check our 
chosen parameters by producing a residual image. While 
we obtain a very good residual image with all prominent 
sources subtracted well, the source detection algorithm 
does not find a source at the target position even for very 
low detection thresholds {2a). Therefore we apply aper- 
ture photometry at the position of the target as derived 

* In principle, we might have detected an as teroid, but the high 
latitude and unusual SED {Trilling ct al."2008') argue against this. 
^ http://ssc.spitzer.caltech.edu/niips/rnipsinstrumcnthandbook 




Fig. 2. — MIPS 24/im mosaic for the position of PSR 
J0437-4715. The image covers 40" X 40". The VLBI radio position 
of PSR J0437-4715 obtained by DcUer et al. (2008) and corrected 
for the proper motion to the epoch of the Spitzer observations is 
plotted as a black cross. The used aperture (r = 3" 5) is also shown, 
together with illustrative overlay contours of the IRAC channel 2 
mosaic. The small white contour south-east of the target marks 
the position where the faint neighbouring source is seen in IRAC 
1 and 2. 

from IRAC channel 1. The MIPS 24 /j,m astrometry gen- 
erally agrees well with the IRAC astrometry and sources 
of IRAC channel 2 are over-plotted in Figure [2] to illus- 
trate this. For the aperture photometry we use a small 
source aperture (3'.'5) and sky annulus (6" — 8") to avoid 
close-by bright neighbors. We note that the aperture 
correction itself has an uncertainty of around 5%, ac- 
cording to the MIPS instrument handbook. The 'mode' 
algorithm for the background estimation in the annulus 
is employed. We derive an aperture-corrected flux den- 
sity of 38 ± 22 /iJy (Icr confidence), or 3cr upper limit of 
105 ^i3y. 

2.2. Magellan NIR photometry 

We observed PSR J0437-4715 with Pe rsson's Auxil- 
iary Nasmyth Infrared Camera (PANIC; iMartini et al.l 
I2004D on 2005 October 20 with the 6.5-m Magellan I 
(Baade) telescope. The conditions were close to photo- 
metric, with seeing of about 0'.'5. We used exposure times 
of 18 X 30s (X^-band), 9 x 20s (H-hand), and 9 x 30s 
(J-band), taken with small pointing dithers between the 
exposures. Data reduction followed standard procedures: 
we masked bad pixels, subracted dark frames and fiat- 
fielded using well-exposed twilight flats. Sky subtrac- 
tion was done using a median-combined image of the 
data. We then combined the individual exposures to- 
gether, identified the visible stars, and masked them for 
a second round of sky subtraction. We then made a final 
mosaic of the individual exposures, see Figure [3] 

We calibrated the photometry relative to 2MASS. 
However, there was not a sufficient number of unsatu- 
rated stars in the field of PSR J0437-4715 to do this 
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in one step. We therefore used other fields observed on 
the same night over a range of airmasses, generahy below 
1.2. Using only stars that were still in the linear regime 
of the detector and were not crowded (1-6 stars in each 
field over 5 fields), we used aperture photometry in IRAF 
with a large aperture of 1'.'9 radius and a sky annulus 
from 2'.'5 to 3'.'8. We corrected to airmass of 1.0 using 
approximate airmass coefficients for the site (0.10, 0.03, 
and 0.07 magairmass^^ for J, and Ks respectively) 
and determined the zero-points for the night. The final 
photometry for PSR J0437-4715 is J = 18.59 ± 0.02, 
H = 18.37 ± 0.06, and = 18.38 ± 0.06, where the 
errors are statistical only. We estimate that the uncer- 
tainties in the photometric zero points were at most 0.03 
mag for Kg and iJ, and 0.02 for J. 

2.3. Optical Spectroscopy 

Spectra of the optical counterpart to PSR J0437— 4715 
were taken on the night of 1999 July 15 using FORSl 
at the Very Large Telescope. The setup wa s the s ame 
as that used by Ivan Kerkwiik fc Kulkarnil (|200l[ ) for 
a spectroscopic study of the faint neutron star, RX 
J1856.5-3754 (Grism 300V with R «300, 2.6 A/pix dis- 
persion; 1" wide long-slit, 0.2"/pix spatial scale). The re- 
duction followed that described by those authors. Briefiy, 
two 45-minute spectra were taken, covering the range of 
3600 to 9000 A at a resolution of ~13A, with the slit po- 
sition an gle chos e n such that both the counterpart and 
star 1 of iBailvnl (|1993l ) were coverec0. The reduction 
involved bias subtraction, sky subtraction, and optimal 
extraction (Hornc 1986) of the spectra of the counterpart 
using the spatial profile of star 1. 

Flux calibration was done in two steps. First, we cal- 
ibrated the fluxes of star 1 using a shorter, 5-minute 
spectrum taken through a wide slit (with the instru- 
me ntal response determined from two fiux standards; 
see lvan Kerkwijk fc Kulkarnill200il especially their Sect. 
4.4). Next, we estimated the slit losses by fitting a 
quadratic function to the ratio of the narrow-slit to wide- 
slit spectra of star 1, and used this to calibrate the 
fluxes of the counterpart. Overall, we believe our rela- 
tive fluxes should be accurate to about 5% between 4500 
and 8800 A. Shortward of 4500 A, the calibration is more 
uncertain, since it is less clear whether our slit losses are 
corrected well: the ratio of the narrow- to wide-slit spec- 
tra of star 1 shows systematic, ~ 10% deviations, likely 
because the wide-slit spectra were taken at high airmass. 
The absolute flux calibration is also less certain, since 
there was some cirrus at the start of the night. We flnd 
that we have to scale the fluxes up by 8% to match the 
HST photometry in §2.4. Furthermore, some regions of 
the spectrum are contaminated by atmospheric absorp- 
tion and high background - we have simply excluded 
these areas from consideration. 

2.4. HST 

We analysed several observations of PSR J0437— 4715 
with the HST, both imaging (WFPC2 and ACS) and 
spectroscopy (ACS). The Advanced Camera for Surveys 
(ACS) observations consisted originally of five visits, one 
in the NUV and optical using the High-resolution Cam- 

^ The star is also known as 086-010575 in UCAC3. 



era (HRC) and the remainder in the FUV, using the 
Solar-blind Channel (SBC) detector. 

2.4.1. WFPC2 Imaging 

We retrieved WFPC2 observations of the field of PSR 
J0437-4715, taken in the two wide-band fihers F555W 
and F814W (roughly V and I band) on 1996-05-19. There 
were four exposures of 140 s each, for a total of 560 s 
taken through each filter. The source was located on the 
PC chip of the array in every case, with a pixel scale of 
0.0455"/pix. P hotometry was performed using HSTphot 
()Dolphinll2000[ ) . a software package optimized to WFPC2 
photometry by profile fitting. It efficiently rejects bad 
data (hot pixels and cosmic ray hits) and fits stellar pro- 
files from libraries of pre-built PSF images (originally 
created using TinyTim). 

For the source, we obtained (Vega-magnitudes) 
?7T-F555W = 20.899(10) and mpsMW = 19.416(8). Since 
the observations were taken early in the WFPC2 lifetime, 
the charge transfer efficiency (CTE) corrections are very 
small. 

2.4.2. ACS Imaging 

Imaging observations are required by the prism spec- 
troscopy, as a way to set the wavelength scale for each ex- 
tracted spectrum. These images are useful, in their own 
right, for accurate broad-band photometry. From our 
observations, we have three bands, F555W, F330W and 
F140LP, each with a detection of the PSR J0437-4715 
system. One exposure was acquired in each single orbit, 
at the start of each orbit for the NUV and the end for the 
FUV. We performed aperture photometry on these, with 
an aperture size set to contain about 90% of the fiux, 
which maximizes the signal-to-noise. The correction to 
total flux was done using the encircled e nergy tables in 
the A CS Instrument Handbook (see also iSirianni et ahl 
l2005f l. and converted to flux using the information sup- 
plied in the image headers. These are based on the latest 
calibration at the time of download (January 2010). Fig- 
ure [3] shows the fields centered on PSR J0437— 4715 in 
the three broad-band filters. 

For F555W and F140LP, multiple images could be used 
to effectively remove cosmic rays and detector blemishes. 
We only had one F330W image to work with, however, 
so such rejection was not possible. No CR is seen in or 
around the photometry aperture. 

The CTE correction for the flux in F555W is very small 
because the source is comparatively bright in that filtei0. 
On the other hand, the F330W fiux does have a substan- 
tial correction of the order 5%-10%. In this case, the 
flux is already relatively uncertain (see Table 3), and the 
CTE effect also increases the uncertainty. 

2.4.3. NUV Spectroscopy 

For the reduction of all the spectral data, we used the 
slitless spectroscopy package aX(fl, Version 2.0 (released 
November 2009), in conjunction with IRAF/PyRAF. A 
detailed description of t he alg orithms used by aXe are 
given in iKiimmel et al.l (j2009l ) . In summary, a spectral 

http: //www, s tscl ■ edu/hst/acs/documents/handbooks/l 

cyclel8/cover .html 

^ www. stsci . edu/hst/acs/documents/isrs/isr0901 .pdf 

^ http: / /www . stecf . org/sof tware/slitless_sof tware/axe/| 
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Fig. 3. — Broad-band imaging of PSR J0437— 4715; north is up. In each case, the source is in the center of the image. Note the cosmic 
rays affecting F330W (because this was a single exposure, see text; we have masked cosmic rays). No counterpart of the faint IRAC 
neighbor of PSR J0437— 4715 (Figure 1) is seen in the ifg-band image. The field of the ACS images is shown on the F814W image. 



TABLE 2 

Broad-band photometry of PSR J0437— 4715. 



Origin 


Band 


Central 
Wavelength^' (A) 


Bandwidth 

(A) 


Spectral Flux F\ 
(10^^* ergs^^cm^^^-i) 


HST 


F555W 


5356 


357 


15.58(15) 




F330W 


3363 


174 


1.9(2) 




F140LP 


1451 


127 


2.93(9) 


Magellan 


J 


12200 


2130 


11.6(4) 




H 


16300 


3070 


5.2(4) 




Ks 


21900 


3900 


1.77(12) 


Spitzer 


IRACl 


35500 


7500 


0.37(6) 




IRAC2 


44900 


10150 


0.17(5) 




IRAC3 


57300 


14250 


0.07(3) 




IRAC4 


78700 


29050 


<0.22 




MIPS24 


240000 


91300 


0.38(20); <1.05 


Danziger et al. 


B 


4420 


940 


8.17(16) 




V 


5450 


880 


16.8(10) 




R 


6380 


1380 


21.1(6) 




I 


7970 


1490 


17.9(7) 



Note. — 



Numbers in parentheses indicate 1-sigma statistical uncertainties in the last digit (s). 
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observation requires a direct (broad-band) image taken 
with the same pointing, which serves to cahbrate the 
wavelength scale on the subsequent dispersed image. For 
the tilted spectrum, pixels within an adjustable distance 
in the orthogonal direction to the trace are summed. The 
background is determined in a region at large distance 
either side of the trace. 

The derivation of the wavelength solution, sensitivity 
curve and aperture correction are detailed in STScI In- 
strument Science Report (ISR) 060^3- 

The ACS HRC is sensitive in the wavelength range 
1700-lOOOOA. For the PR200L prism observations, we 
used the direct images taken in the F555W filter. Since 
all the exposures for the NUV were taken during a single 
visit, they are very well aligned with one-another. Al- 
though in slitless spectroscopy spectra can in principle 
overlap, in our case, the source is well-isolated and the 
brightest in the field, so such contamination is not an 
issue. 

Although the individual exposures are the ones used in 
the analysis, running multidrizzle ensures the masking 
of cosmic rays as bad data, and provides an opportunity 
to perform visual inspection of the data. Cosmic rays 
are not removed, but a fiag is set in the data quality 
(DQ) extension of the input image, marking that pixel 
as unusable. In order to make use of this filtering out 
of bad data, it is necessary to copy the modified DQ 
extension into the original, unmodified image file, and 
use this as input for aXe. 

Figure |4] shows the combined spectral image in the 
vicinity of PSR J0437— 4715. A number of interesting 
features are seen. The bright elongated source is the 
area where the dispersion becomes very small, and all 
the light for wavelengths A > 5000 A ends up centred 
on just a few pixels (the red pile-up effect). Since the 
source is intrinsically red (the cool WD dominates, see 
§3.1.2), this large halo dominates over the faint spec- 
tral trace. Faint discrete structures are seen around the 
bright source, probably associated with reflections within 
the optical system, and restricted to A > 4000 A. A dif- 
fuse halo of greater extent is also seen, possibly the wings 
of the PSF in the cross-dispersion direction. Not seen by 
eye, but present in the background analysis, the bright 
core at pixel (800,570) also has a diffraction spike along 
the X-axis, at a small angle to the spectral trace. This 
limits the size of extraction region which will be free of 
conta mination, particu larly for intrinsically red sources 
(e.g., iWalkowicz eTall llJOSV 

We calculated final fluxes by heavily binning the count 
rates measured by aXe, together with the provided sen- 
sitivity curve and aperture correction (see ISR 060^0) • 

We find, after much experimentation, that the op- 
timum trade-off between gathering enough signal and 
gathering too much noise and contamination was at an 
aperture half- width of IxFWHM (7pix fuU width). For 
the background estimate we employed two regimes: at 
the faint, short- wavelength end (A < 3175 A), where 
there was no bright halo, and the diffraction spike would 
be far from the source, we used a region from 20xFWHM 
to 21.5xFWIIM from the trace. For longer wavelengths. 



" see http://www.stecf.org/documents/isr/isr0603.pdf 
see http://www.stecf.org/documents/isr/isr0603.pdf 



where the halo contamination will be significant, we used 
a background sampled from 1.5xFWIIM to BxFWHM 
from the trace in order to subtract this contribution, but 
corrected the fiux upwards for the small amount of source 
flux («6%) that will fall in the background sampling re- 
gion. We also compared the results with those obtained 
using optimal extraction. In theory, optimal extraction 
should enhance the signal-to-noise ratio if the trace is 
well sampled in the spacial direction. We found that 
the optimal extraction produced a smoother spectrum, 
in better agreement with the broad-band photometric 
point at 3363 A, so we have reason to believe that it is a 
fair extraction in this case; we therefore use the optimal 
extraction fluxes. 

The final uncertainties on the fiuxes are a combination 
of Poissonian photon noise (in both the signal and the 
background), small uncertainties in the sensitivity curve, 
and larger systematic calibration uncertainties of 5% in 
the range 1800-3000 A rising to 15%-20% at the edge of 
the sensitivity range (Kiimmel, pers. comm.). All the 
uncertainties were added in quadrature. The final set 
of counts and fluxes for the wavelength bins chosen are 
listed in Tabled 

2.4.4. FUV Spectroscopy 

PSR J0437-4715 was imaged through the PR130L 
prism. This prism was chosen in order to decrease the 
background level, since it has low throughput shortward 
of 130 nm, and thus limits the contribution of the bright- 
est geocoronal lines. Furthermore, the spectral expo- 
sures were taken while the spacecraft was in the Earth's 
shadow, further reducing the background. The imag- 
ing exposures were done in the remainder of the orbit. 
In addition, the SBC's Multi- Anode MicroChannel Array 
(MAMA) detector does not suffer from cosmic ray hits 
as do CCD detectors. As the resulting combined image 
in Figure [S] shows, the background is faint and uniform, 
and the spectrum is very easy to see. In comparison to 
the NUV observations, the spectral trace is well-deflned, 
and there is no background structure or any halo/pile- up 
effects. 

Reduction proceeded in a similar way to the NUV 
spectral analysis above. Since the exposures come from 
separate visits in this case, and therefore required re- 
acquisition of a guide star, alignment between them is 
not guaranteed, and the source position on the direct 
image was translated back to each input direct image 
using the aXe task lOLPREP. This was particularly im- 
portant for the Visit 02 exposures, because of the loss 
of the guide star during this visit. The data had an 
enhanced background level, and the data from the last 
orbit were totally lost. Nevertheless, we do include the 
data from the first part of the visit, and we find that it 
improves rather than decreases the signal-to-noise ratio. 
Before proceeding with the extraction, the input coor- 
dinates need to be corrected for an offset caused by the 
choice of broad-band fllter for imaging (see ISR 060Q. 

Due to the much lower background level compared to 
the NUV, and absence of the problems associated with 
red pile-up and diffraction spikes near the spectral trace, 
it proved to be optimal to extract the source from a some- 
what wider aperture, half-width of 4xFWHM (19 pix full 

|http: //www. stecf . org/documents/isr/isr0602 .pdf | 
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Fig. 4. — NUV spectral image of PSR J0437— 4715 through prism PR200L. Wavelength decreases left-to-right, non-linearly, as marked 
(values in A). The large source seen at long wavelengths (left) is due to the red pile-up, see text. The upper panel shows the bright red 
source caused by the red pile up, and its surrounding reflection and diffuse structures. The lower panel is the same, but with the color scale 
stretched, and the extraction regions marked. The background regions used for long wavelengths are also marked; the background regions 
for short wavelengths lie far from the spectrum (sec text). 



width). See Figure [S] for the locations of the extraction 
region. The results are listed in Table [31 With the low 
sensitivity, the systematic uncertainties in the calibra- 
tion make less of an impact in this case, and optimal 
extraction produced a spectrum that looks a little nois- 
ier, but otherwise similar - we therefore used simple fixed 
aperture extraction. Systematic uncertainties are of the 
order 5% throughout the spectral band analyzed. 

2.5. X-rays 

We retrieved one long XMM-Newton observation and 
nine Chandra observations from the archives, and re- 
processed all the data to be analyzed simultaneously. 
The data consist of CCD imaging and provide low- 
resolution spectra using the inherent energy resolution 
of the detector, as opposed to inserting a spectral ele- 
ment such as a grating. Although the Reflection Grating 
Spectrometer (RGS) instrument of XMM-Newton does 
provide such data, its count rate was too low for consid- 
eration. 

2.5.1. XMM-Newton 

XMM-Newton observed PSR J0437-4715 on 2002-10- 
09, 03:32:34-22:50:04. This observation and its analysis 
are described in detail by Zavlin (2006). Here we consider 
the data from the two EPIC-MOS detectors, operated in 
Full Window mode (2.6 s time resolution) and the EPIC- 
PN detector operated in Fast Timing mode (0.03 ms time 
resolution). The raw data were processed into calibrated 
event files using the standard tasks of the provided soft- 



ware XMM-SAS, version 9.0 (20090615_1801). After fil- 
tering the data for times of flaring, the total live time 
was 63.0 ks for the MOSs and 62.9 ks for the PN. 

For the MOS data sets, we extracted the source spec- 
trum from a circle of 15" radius centered on the source. 
This only captures some 70% of the total flux, but such 
a compromise is necessary because of the relatively high 
background and the background from a circular annulus 
with inner and outer radii of 114"and 200", respectively. 

In timing mode, all information about the Y-axis po- 
sition of a photon on the detector is lost. We filtered 
out low-energy electronic noise[3 by keeping only ener- 
gies E > 600 eV. There is a bright AGN in the field whose 
one-dimens ional PSF partly overlaps with that of the 
pulsar (see iZavlin et al.l I2002D . We therefore extracted 
counts from an aperture chosen to avoid contamination, 
that only contained 61%-66% of the flux (higher frac- 
tion for lower energies). We correct the spectrum for 
this finite-aperture effect. 

2.5.2. Chandra 

We retrieved nine Chandra ACIS imaging datasets, 
with IDs 6154, 6157, 6767, 6768, 741, for a total ex- 
posure time of 73.2 ks. We do not consider High Resolu- 
tion Camera (HRC) observations nor ACIS observations 
affected by pile-up. The observations have a range of ex- 
posure times and off-axis distances. The total exposure 
time was 73.2 ks. 

1^ http : //www . star . le . ac.uk/~amr30/BG/BGTable .html ", col- 
umn titled Electronic Noise 
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TABLE 3 

FUV AND NUV SPECTRAL DATA 



Wavelength bin 


Bin width 


Gross Counts 


Background Counts 


Net Counts'' 


Flux^ Fx 


centre (A) 


(A) 








(10~^* ergs~^cm~^A~^) 


1245 


30 


403 


225 


178(20) 


6.42(0.80) 


1275 


30 


424 


177 


247(21) 


5.34(0.53) 


1305 


30 


412 


144 


268(20) 


5.19(0.48) 


1335 


30 


387 


106 


281(20) 


5.36(0.47) 


1365 


30 


348 


93 


255(19) 


5.00(0.45) 


1395 


30 


256 


70 


186(16) 


3.86(0.39) 


1425 


30 


253 


70 


183(16) 


4.05(0.41) 


1455 


30 


171 


51 


120(13) 


2 86f0 35^ 


1485 


30 


188 


49 


139(14) 


3 fi0?0 40^ 


1515 


30 


195 


49 


146(14) 


4 1 fi?0 4fi^ 


1545 


30 


132 


39 


93(11) 


2 Qq/o 40^ 


1575 


30 


128 


41 






1605 


30 


119 


39 


80(11) 




1645 


50 


163 


55 




3 27(0 4'Vi 


1695 


50 


107 


41 


66(10) 


2 42f0 41 \ 


1745 


50 


111 


41 


70(12) 


3 36 fO 55^ 


1795 


50 


99 


41 


58(10) 


4.26(0.78) 


1845 


50 


53 


27 


25(7) 


3.16(0.94) 


2150 


700 


54808 


52685 


2123(230) 


1.54(0.64) 


2600 


200 


7699 


7006 


693(88) 


1.07(0.15) 


2750 


100 


3053 


2648 


405(55) 


1.27(0.18) 


2850 


100 


3127 


2731 


396(56) 


1.15(0.17) 


2950 


100 


3167 


2516 


651(56) 


1.72(0.17) 


3050 


100 


2177 


1689 


488(47) 


1.17(0.13) 


3150 


100 


2620 


1894 


726(51) 


1.68(0.15) 


3300 


200 


4635 


3243 


1392(68) 


1.73(0.27) 


3500 


200 


4096 


2377 


1719(64) 


2.33(0.36) 


3700 


200 


2506 


1203 


1303(50) 


1.92(0.30) 


3900 


200 


2824 


1238 


1586(50) 


2.54(0.41) 



" Uncertainties are statistical only for the net counts, and include both statistical and systematic/calibration contributions for the flux. 
Numbers in parentheses show one-sigma uncertainties. 
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Fig. 5. — FUV combined spectral image of PSR J0437— 4715 through prism PR130L, as produced by multidrizzle. The extraction 
region is labeled, along with approximate wavelengths in A. 
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The source extraction region was chosen to be fairly 
large and uniform for all the observations, as the back- 
ground is low, and the varying ofF-axis distance for each 
observation means the size of the PSF is different. The 
radius was set at 9.15 pix (4'.'5), which includes some 99% 
of the flux for the on-axis PSF, and most (w90%) even 
for the broadest PSF in the set of observation. The back- 
ground was taken either from an annular region centered 
on the source, or rectangular regions either side of the 
readout direction, for observations with subarrays, where 
the active area is restricted in one direction. 

2.6. Fermi 

We retrieved data on the region containing PSR 
J0437-4715 collected in 18 months of Fermi LAT ob- 
servation. 7-ray data has much lower count-rates than 
X-rays, extended, broad-winged PSF and frequency- 
dependent energ y resolution. Standard filtering (see 
lAbdo et"ani2010D was applied. The reduction to obtain a 
spectral model or fluxes for a single source is therefore a 
likelihood maximization method, with explicit considera- 
tion of nearby sources and their respective spectra. Since 
spectral fitting is then an intrinsic part of the extraction 
process, we describe it in Section [5751 

3. SPECTRUM AND FITS 

In Figure [S] we show the UV-IR spectrum of J0437. 
Note that we don't include the MIPS point, since it is 
more likely to be a (non-restrictive) upper limit. It is re- 
assuring to see the high degree of agreement between our 
FUV spectral points, the F130LP broad-band photome- 
try and the previous STIS spectrum. Likewise, the NUV 
spectrum and the F330W photometry agree, but there 
is no independent data in this region to verify the spec- 
tral shape. Furthermore, the NUV and optical points 
support a smooth continuum linking the two spectral re- 
gions. Given the several problems inherent in processing 
the PR2G0L data, we regard with suspicion the hint of 
structure around 3700 A. 

The extinction to the source is not well constrained by 
fits to the X-ray spectrum, but must be s mall, with the 
hydro gen column of ^ 10^° cm~^ or less (jZavlin et al.l 
I2002D . This is equivalent to E{B - V) < 0.05, consis- 
tent with the small distance to the source, d = 156.3 pc. 
We therefore attempt to model our combined spectrum 
with extinction values E{B - V) = 0.01, 0.03 and 0.05. 
Furthermore, there is some considerable difference in the 
extinction curves in the literature, for a given value of 
extinction. This is an additional uncertainty in our fits. 
In practi ce, the shapes of th e extinction cu r ves in , for 
instance, iSeatonI (|1979f ) and iCardelh et"all ()1989[ ) are 
very similar, but diff erent in normalization. We will use 
iCardelli et al.l (|1989| )'s more recent parameterization for 
our fitting. 

In the NUV and FUV we considered only the spectro- 
scopic fluxes because of the lower signal-to-noise of the 
F330W point and the extreme asymmetry of the F130LP 
band-pass. 

3.1. UV to mid-IR 
3.1.1. Analytic models 

We first fit the observed spectrum with simple ana- 
lytical models. In particular, the maximum in the op- 
tical can be fitted with a black-body (BB) curve, and 



TABLE 4 
Fits to IR-UV spectrum 



E{B - V) 


T (K) 




a 


A^pl'^ 






BB-I-PL 






0.01 


3870(40) 


2.02(11) 


-3.3(4) 


3.7(2) 


0.03 


3910(40) 


1.97(11) 


-3.4(4) 


4.3(3) 


0.05 


3950(40) 


1.95(11) 


-3.5(4) 


5.0(3) 






WDA+PL 






0.03 


3950(150) 


1.87(20) 


-2.46(17) 


4.61(19) 



Note. — Black-body plus power-law fits (§3.1.1) and WD atmo- 
sphere model plus power-law fit (§3.1.2). Uncertainties are 1-sigma, 
taken from the diagonal of the covariance matrix and assuming the 
fit is reasonable. 

^ Effective radius for a distance of 156.3 pc, units of 10* km. 

^ Spectral flux of the PL component at 1500 A, in 

10~^® ergs^^cm^^A^l. 

the FUV-NUV with a power-law (Fa oc A"; PL). We at- 
tempt this for different reasonable values of reddening, 
E{B - V) =0.01, 0.03, 0.05. The results of these fits 
are shown in Table 21 It should be noted that none of 
these fits is statistically acceptable (see Figure [7]), but 
the numbers can be taken as indicative. The spectrum 
shortward of 1500 A appears to be well-represented by 
the a K, —3.5 power-law. 

We can see that a simplistic sum of a ~4000K black- 
body and a PL with a — —(3 — 4) (i.e., close to a 
Rayleigh- Jeans tail of a thermal spectrum with a high 
temperature) roughly describes the data. The apparent 
BB emitting radius is about 20,000km (or 0.028 i?©). 

3.1.2. WD atmosphere 

Although informative, the BB only crudely approx- 
imates the optical-IR spectrum above and cannot re- 
veal information on the surface properties of the WD. 
A better description could be provided by a detailed 
white dwarf (WD) atmosphere model. Atmosphere mod- 
els have been developed fo r a variety of WD m asses, 
ages and compositions (e.g.. iBergeron et "aLlll995D . Note 
that, in principle, the WD surface is irradiated by the 
higher-energy emission of the companion pulsar, but 
the effect is small: for the known binary separation 
(jVerbiest et al.l 120081) , the radiative intensity from the 
pulsar (see also §4) at the distance of the WD is approxi- 
mately 3xl0®ergs~^cm~^, about 2% of the outward ra- 
diative fiux; we therefore neglect it in our modeling. We 
describe in detail the model fitting process. 

The spectrum of the WD companion was modeled us- 
ing white dwarf atmosphere (WDA) models that account 
for recently introduced improvements in the description 
of the atmosph ere of these stars, especial ly Ly-a absorp- 
tion in the UV (IKowalski &: S aumon' 2006), op acity of he- 
lium dominated mediu m dKowalski 2006; Kowalski et all 
I2007D and refraction (jKowalski fc SaumonI I2004D . The 
models used in the analysis were successful in describ- 
ing the entire spectral energy distributions of many cool 
WDs, including the UV spectrum of pure-H atmosphere 
WD BPM 4729 (jKowalski k Saumon 2006 ) a nd the near 
UV s pectra of WDs with T^s ~ 4000 (jKihc et al.l 
I20T0I) . 

The fit in effective temperature and gravity for He and 
H atmospheres was performed for the known photometry, 
BVRI of ,Danziger et al.. (1995) . and Magellan JHK and 
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Fig. 6. — UV to mid-IR spectrum of PSR J0437— 4715 with a zoom of the UV portion in the inset. Black points are from prism 
PR130L (UV) and VLT/FORSl (optical), red points from prism PR200L, blue points from HST broad-band imaging, cyan triangles from 
Magellan/PANIC IR photometry, and yellow squares (and a limit) from Spit zer photometry. Cyan UV points are from previous STIS 
spectroscopy (K04), and green points are ground-based photometry from (Da nziger et al.lll993l ). Open markers denote fluxes from FORSl 
or PR200L with uncertain calibration. 



Spitzer IR measurements reported h ere, assuming that 
the distance to the system is 156.3 pc (jPeller et al.ll2008l) 
and the reddening is E{B — V) = 0.03. The spectral en- 
ergy distribution of the WD is best reproduced by a pure 
hydrogen atmosphere model. The hydrogen is clearly vis- 
ible, when comparing with pure-He model spectrum, by 
strong a bsorption in the UV, due to strong Ly-a red wing 
opacity (iKowalski fc Saumonll2006n and flux suppression 
in the IR due to collision-induced absorption (CIA) by 
molecular hydrogen (Saumon & Jacobson 19990. The 
best-fit model closely matches the VLT/FORSl spec- 
trum (see Figure [7]). 

The best-fit effective temperature, T = 3950 ± 150 K 
is in agreement wi th previo us estimates (iBell et al.lll993l : 
iBassa et all 120061 DanzigeretaD |1993D and Ihe BE 
fit from above, see Figure [7l The obtained gravity 
(logg = 6.98 ± 0.15 [cgs]) fits the models of a helium 
core white dwarf (jSarna et al.ll2001l : IAlthaus et al.ll2009D . 
assuming the mass of the companion to be 0.254 M0 
(jVerbiest et al.l [20081) . The uncertainties are estimated 
from the uncertainty in the photometric fluxes (Table 
[2|), and the small uncertainties in the distance and WD 
mass. We include the fit parameters in Table H] we also 
list the radius, R, which is fitted implicitly via log g (be- 
cause the mass is well known); the quoted error reflects 



uncertainties in the temperature and distance. With- 
out spectral lines, the spectrum provides no distance- 
independent constraints on the gravity). 

While some of cool WD pulsar companions have optical 
photometry, this is usually not enough for definitive as- 
signment of their atmospheric chemical composi tion and 
parameters (see discussion in IBassa et al.|[2006D . At so 
low temperatures the admixture of helium should pro- 
duce substantially different CIA effect due to H2-IIe col- 
lisions and therefore would be different from pure hydro- 
gen atm osphere spectrum in the IR, as indicated in Fig- 
ure 5 of iHall et al.l ()2008D . In this study we present the 
full spectral energy distribution of the companion, and 
by obtaining good match with the pure-H atmosphere 
model we can be confident that the object is indeed a 
H- atmosphere WD. This allows for proper derivation of 
the atmosphere parameters (Toff, g) and the age of the 
WD. Measurements of near- and -mid IR fluxes are es- 
sential for a definitive assignment of the atmosphere com- 
position of other cool WD pulsar companions. The WD 
companion atmosphere of PSR J0751-I-18 07 was classi- 
fied as helium-rich by Bas sa et al.l (|2006l ). which is at 
odds with the evolutionary models for such stars. In 
that particular work such a conclusion was reached by 
performing analysis with the WD atmosphere models of 
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Fig. 7. — Top: Black-body plus power-law fit to the spectral data for PSR J0437— 4715. This is the case for the preferred value of 
reddening E{B — V) = 0.03, giving Tbb = 3910 K, a = —3.4. The lower panel shows fractional residuals (flux-model)/model. Bottom: 
Model with white dwarf atmosphere (T(,ff = 3950 K, logg = 6.98, pure H) and power-law, with index a = —2.46 (solid blue). The two 
components are shown separately (dashed blue). The lower panels show fractional residuals; (fiux-model)/model. 
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iBergeron et al.l ()1995[ ). However, these models tend to 
bias the prediction of the atmospheric content of a cool 
white dwarfs from hydrogen-rich to hehum-rich , as they 
do not fully account for the red wing of the Ly-a ab- 
sorption and the free-free opacity of dense helium (see 
discussion in Kowalski & Saumon (2006). Our conclu- 
sion for the atmospheric composition of WD companion 
of J0437-4715, obtained by a good fit to the entire WD 
spectral energy distribution by the pure-H atmosphere 
model, is in line with other works, which postulate a 
pure hydrogen ou termost layer (Hansen & Phinncv 1998; 
ISarna et alJl200l . 

Using this (now fixed) model for the WD spectrum, 
we can again fit a PL to the pulsar part of the spec- 
trum. We choose to fit the range A < 3575 A, where 
the pulsar dominates. The fit is to be compared with 
the black-body fits in Table H we find a = -2.46(17), 
Npi, = 4.61(19) X IQ-i^ergs-^cm-^A-^ From Figure 
[7l it would appear that the fit is good, except for the 
portion of the spectrum in the NUV, where the WD at- 
mosphere and the PL are of similar strengths. Also, the 
value of the power-law index is somewhat puzzling, since 
it is no longer thermal-like, nor is it typical for a pul- 
sar non-thermal spectrum. Indeed, the residuals in the 
FUV appear to show a systematic slope. This suggests 
that the UV spectrum may actually consist of two com- 
ponents, for instance a thermal Ray leigh- Jeans tail (PL 
with a = —4) and a non-thermal PL more typical of 
old MSPs (PL with -1 > a > -2, i.e., photon index 
1 < r < 2) - see Section 4. 

3.1.3. Thermal FUV 

Using HST/STIS, K04 were able to measure the FUV 
spectrum in 4 bins with low S/N, which rendered the 
measurement of the spectral slope rather uncertain. Our 
much better quality FUV spectrum from ACS matches 
those STIS fluxes rather well (Figure [S]) and also the 
FUV broad-band photometry. This gives us confidence 
in the accuracy of the cross-calibration between the in- 
struments. 

Under the assumption that the FUV PL-like emission 
is indeed a Rayleigh- Jeans law {F\ oc A^^ or Fi, oc v'^), 
we can constrain the temperature and radius of the emit- 
ting body for the known distance, d = 156.3 pc. We cal- 
culate Ti?^ ^ 1.9x lO'^Kkm^. In the case here, however, 
a black-body law starts to deviate from a Rayleigh- Jeans 
law in the FUV for temperatures T < 5 x 10^ K. Note 
that R = i?oo and T = Too, i-e., the temperature and 
radius as measured by a distant observer. 

We show possible thermal fits to the FUV spectrum 
in Figure [H We include the lowest-energy points from 
the X-ray spectra in Section 13.21 - these constrain the 
high-temperature limit of the possible thermal emission 
toT < 3.5xlO^K. The lower limit on the temperature is 
set by requiring that the emitting area be reasonable for 
a NS. Whereas temperatures as low as 7 x 10'' K can still 
roughly fit the data, the required radius, 24 km, is larger 
than the largest plausible value of 15 km (Lattimcr 2010). 
For R — 15 km, the corresponding surface temperature 
is T = 1.25 X 10^ K, and the likely NS size of i? = 13 km 
yields a surface temperature of 1.5 x 10^ K. 

All of these estimates are affected by the choice of 
reddening in the FUV and extinction in X-rays (for 
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Fig. 8. — Thermal fits to the de- reddened FUV spectrum (points 
and error bars, left), showing the coolest temperature allowed by 
the data (T = 7 X 10* K) but with a radius too large for a NS; the 
coolest temperature compatible with a NS-sized emitter ol R = 
15 km, T = 1.25 X 10^ K; and the hottest temperature allowed by 
the X-ray data (extinction-corrected points and error bars, right), 
T = 3.5 X 10^ K with an emitting radius of 7.8 km. These values 
correspond to UV reddening E{B — V) = 0.03 and X-ray extincting 
column Nh = 7 x 10^^ cm"^. 

the upper limit). All the temperatures quoted are for 
E{B — V) — 0.03, and should be corrected upwards 
or downwards by about 15% for E{B ~ V) = 0.05 and 
E{B - U) = 0.01, respectively Similar to the WDA+PL 
model in the previous section, any WDA-|-thermal model 
leaves noticeable residuals in the NUV range (see Figure 

3.1.4. Emission Line? 

In Figure [7] enhanced flux seen in the residuals around 
1800 A, attributable to three bins of the FUV spectrum. 
Figure [5] shows a zoom of the combined spectrum im- 
age, with enhanced contrast, and the vicinity of the line 
marked. No feature is seen in the raw data by eye; the 
spectral response is, however, rapidly falling to longer 
wavelengths in this region. If the point at 1795 A were 
omitted, there would be no question of a line; as it is, 
we cannot rule out for certain the existence of an emis- 
sion line here. We do not, however, attempt to fit a line 
to the three points, just note that their presence may 
slightly affect our accuracy in measuring the slope of the 
spectrum 1000-1700 A. 

No evidence is found of a spectral line at 1372 A, sus- 
pected by K04. Although there does happen to be an 
upward fluctuation in the gross counts spectrum at ap- 
proximately the right wavelength, it is consistent with 
a fluctuation in the background, and smaller in the fi- 
nal flux spectrum than cither other fluctuations or the 
uncertainty expected for Poissonian statistics. As K04 
claimed only a marginal detection, we assume it was due 
to a statistical fluctuation only. 

3.2. X-ray band 

We fitted the set of X-ray observations simultaneously 
using the X-ray analysis package Sherpa. The eight sepa- 
rate spectra were loaded, together with their correspond- 
ing response and background data. Each spectrum was 
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Fig. 9. — Image of the PR130L spectrum, in the region of the possible line at 1800 A. The approximate extent of the feature is marked 
by white dotted lines. 



grouped into bins with a signal-to-noise ratio of least 8 for 
XMM and 5 for Chandra. All fits were performed on the 
whole data-set simultaneously, with energy bounds 0.1- 
10 keV for XMM MOS, 0.6-10 keV for XMM-Newton PN 
and 0.3-8 keV for Chandra ACIS data. 

The results of fits with a variety of continuum mod- 
els is shown in Table El Among others, we used 
the neutron star atmosphere model implemented in 
Xspec as "N SAgrav" for non-magnetised atmospheres 
( Zavlin et al. 1996). Keeping the mass fixed at 1.76 Mq 
( Verbiest et al. .2008i) , and the radius at the default 



10 km, we find a marginally better fit than with a sin- 
gle BB. The hot component best-fit normalization cor- 
responds to a large distance (d > 5 kpc) if the emission 
were coming from the entire visible surface. Since the 
true distance is d = 156 pc, this suggests that the emis- 
sion is coming from a region with effective radius less 
than 3% of the neutron star radius. 

The upper limit established by iZavlin fc Pavlovl 
(|1998f l on the column density to PSR J0437-4715 is 
1.2x 10^° cm~^, based on the X-ray spectrum of the 
nearby AGN. This would correspond to E(B -V) « 0.02 
optical reddening (jPredehl fc Schmit3 119951 ) (note that 
this relation shows considerable scatter, especially for 
low extinction). From fitting the softer EOSAT/PSPC 
spectra with black-body, power-law or neutron star at- 
mosphere models yielded column densities in the range 
(1-9) X 10^^ cm~^. We performed two sets of fits, one 
with the fixed hydrogen column density of Nh — 7 x 
lO^^cm"^, and a second set with the extinction free. 
We u sed the cross-sections of iMorrison fc McCammonI 
(|1983'), as implemented in the model wabs. The effect of 
extinction for E > 0.5 keV is very small. The fit param- 
eters we obtain are generally similar between the fixed 



and free extinction parameter. In Figure [TO] we show 
one of our fits, with the NSAgrav-|-PL model and fitted 
extinction. 

In Figure [TT] we show the average fluxed spectrum (see 
Section 4), plotted with the best-fit BB+PL model. We 
can see that generally the PL component's photon index 
is not only determined by the high-energy tail, but by a 
combination of this tail and the low-energy, extinction- 
affected data. Since the PL cannot rise so steeply into the 
EUV/FUV, we suspect that the X-ray spectrum is better 
modeled by multiple thermal components and a fiatter 
PL. Furthermore, although the fit is formally acceptable, 
the majority of the is contributed by the high signal- 
to-noise 0.5-2 keV range, leaving some structure in the 
residuals. 

The variety of models tested in Table [51 all produce 
statistically acceptable fits, except BB+BB. We see ther- 
mal emission dominating the high signal-to-noise bulk of 
the spectrum with possibly further cool thermal temper- 
atures at the low energies and PL at high energies. If 
we perform the fitting only for the higher energies, we 
find increasingly small (but poorly constrained) values 
of the photon index; e.g., for BB-I-PL fitted above 1 keV, 
r — 2.2(3), and for PL only fitted to energies above 3, 
4, and 5keV, we find F = 1.6(4), 2.0(6) and 1.6(9), re- 
spectively. Since the two-temperature atmosphere mod- 
els also yield similar values of the photon index (Table 
E]), it is more likely a more reallistic description of the 
high-energy spectrum. 

3.3. Gamma-rays 

We modeled the spectrum in the Fermi band by direct 
maximum-likelihood analysis of the calibrated photons 
falling near the position of PSR J0437— 4715 on the sky. 
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Fig. 10. — Count-rate spectrum of PSR J0437-4715, with NSAgrav+PL model and Nh = 7.4 X lO^^cm'^ extinction, which had the 
lowest reduced chi-squared value of our fits. The black points are PN data, green is ACIS and red is MOSl. 




Fig. 11. — The X-ray BB-(-PL model, compared to the combined fluxed spectrum (i.e., each flux point is the weighted average of the 
fluxes from the three XMM-Newton and five Chandra spectra in that energy bin). The model was obtained by standard 'forward fitting' 
of the count rate spectra. Fractional residuals are shown in the bottom panel (i.e., (f-model)/model). 
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TABLE 5 
X-RAY Spectral Fits 



Model 


Nh 
(IQlS cm-2) 


kTi 
(eV) 


(km) 


kT2 
(eV) 


R2 
(km) 


r 




xVdof 


BB+PL 


7 


216(2) 


0.073(3) 






2.699(19) 


1.28(3) 


646/629 


BB+BB 


7 


286(5) 


0.045(5) 


108(3) 


0.35(4) 






1084/629 


BB+BB+PL 


7 


258(12) 


0.044(7) 


141(13) 


0.112(15) 


2.69(4) 


1.05(4) 


598/627 


NSAgrav+PL 


7 


138(2) 


0.294(12) 






2.69(3) 


1.05(3) 


601/629 


NSAgrav+NSAgrav+PL 


7 


140(20) 


0.4(2) 


30(11) 


7(2) 


1.9(6) 


0.31(2)l«-6 


599/627 


BB+PL 
BB+BB+PL 
NSAgrav+PL 
NSAgrav-NSAgrav+l'L 


13.9(13) 

7(3) 


223(3) 

24ot5E; 

110(14) 
l.-)8(13j 


0.065(2) 

'J 0.04 

0..36(2) 
0.2 !(.-,) 


125+f 

18(9) 


0.15+°;f 

2,2(8) 


2.86(3) 
2 6q+"''^ 

^■"'^-0.08 

2.68(5) 
1,()(1) 


1.46(4) 
1.1(2) 
1.05(9) 
0,23(19) 


613/628 
599/626 
•599/628 
592/()26 



Note. — PL: power-law, BE: black-body, NSAgrav: neutron star hydrogen atmosphere (non-magnetized). Numbers in parentheses 
show the one-sigma confidence interval in the last digit (s). Models in the upper section have fixed column density, Nh, and in the lower 
section it is allowed to vary. For NSAgrav, T = T^ffy the model effective temperature (uncorrected for red-shift) and the radius R of the 
effective emitting area. 

^ PL normalization A in units of 10~* phs~^cm~'^keV~^ at IkeV. 



17 



In order to derive spectral models and fluxes for 
J0437, we followed the standard procedure for Fermi 
LAT datgp^. In particular we modeled out emission from 
IFGL sources (see Abdo et al. 2010a) within a r = 10° 
circle centered on J0437. We derived the spectral fit to 
the J0437 spectrum in 0.1-4 GeV ranged by maximiz- 
ing the unbinned likelihood with the gtlike procedure. 
Standard background models were used to account for 
both Galactic and extragalactic diffuse emission as well 
as instrumental background. We have fitted both PL and 
cutoff PL [i?^'" exp(— £'/£'cut)] models. Due to the lim- 
ited statistics available in the data, the cutoff energy was 
not well constrained by the fitting, with a best-fit value 
of iJ(.ut=l-3 GeV, but with a broad acceptable range. We 
obtain best-fit P = 2.6 ± 0.2 and ~ 2.2, for the PL fit 
and the cut-off PL fit, respectively. The observed fiux in 
0.1-100 GeV is 2.7(8) x 10"" ergs'^cm-^ (compared to 
the IFGL value 2.4(7) x 10"" ergs-^cm-^). With the 
statistics available it is not possible to firmly discrimi- 
nate between the PL and cut-off PL models, but better 
S/N sppectra of other MSPs su ggest that cut-offs above 
1 GeV are ubiquitous in MSPs (|Abdo et al.ll2010[ ). 

The slight differe nce between the fi t parameters here 
and those listed in lAbdo et"al] ()2010f ) are likely due to 
the larger amount of data which we analyzed. 

In order to produce fluxes and uncertainties for the 
multi-wavelength fltting below, we spit the available en- 
ergy range into five bins, and repeated the gtlike pro- 
cedure for each bin, assuming cut-off PL. 

4. MULTI- WAVELENGTH ANALYSIS 

We wish to consider whether the ensemble of data 
across many orders of magnitude in photon energy can 
be fit by simple models. 

We converted each spectrum into spectral energy fiux, 
F^. This involves dividing the number of counts in a bin 
with the effective area at that particular frequency and 
multiplying by the bin energy. In the optical/UV regime, 
the effective area is simply the sensitivity function (also 
called zero point for photometry), but at higher energies, 
one normally also considers spectral redistribution: the 
number of counts in a given energy bin is determined 
by the incident spectral flux not only at that energy but 
also at other energies. Whereas in the optical this is 
typically only from the neighboring bins (and sets the 
sample spacing in spectral space), for X- and 7-rays, a 
photon of given energy may generate signal at a very 
diffe rent energy. We define the flux F in bin h as (see 
e.g.. iNowakll2005h : 



Ffiuxcd{h) 



[C{h) - B{h)] ltd 
J R{h,E)A{E)dE 



(1) 



where C{h) and B{h) are the gross counts and back- 
ground in bin R{h, E) is the redistribution matrix i.e., 
the normalized probability of a photon of energy E be- 
ing detected in bin h, and A{E) is the effective area for 
energy E. 

In order to combine the X-ray spectra, we sample the 
flux from each data-set in a given energy bin, and take 

See |http://fe rmi. gsfc.nasa.gov/ssc/data/analysis/scitools/ 1 
There are virtually no photons detected from J0437 above 4 

GeV. 



the mean weighted by the uncertainties. We note that 
although redistribution undoubtedly affects the appear- 
ance of narrow features, for the wide spectral range and 
broad bins considered here, it should make little differ- 
ence. 

In Figure [12] we show the multi- wavelength spectrum 
of PSR J0437-4715 and the fits that were made for the 
spectral bands separately: the IR-UV, X-rays and 7-rays. 
The fits are combinations of PLs and thermal contribu- 
tions, but we show only the PL components, extrapo- 
lated over many orders of magnitude, because a thermal 
component fitted to the X-rays, for example, does not 
contribute in any other spectral window. For the 7-ray 
fit to the Fermi data, there is significant uncertainty in 
the best-fit parameters. In order to calculate the range 
of spectra consistent with the models, we perform 10000 
simulations of spectra, taking parameters from Gaussian 
distributions with central value and spread given by the 
fitting; dashed lines in Figure [T^] show the best-fit mod- 
els, and dotted lines represent the one-sigma confidence 
bounds. For extrapolating the X-ray and optical spectra, 
the PL photon indices are very well defined for a given fit, 
but fits with different models give different, inconsistent 
values. We plot the best-fit PL components of the vari- 
ous models, P = 3, 1.5 for X-rays and a = —4, —3.4, —2.5 
for the UV. Here the value a = —4 corresponds to a 
Rayleigh- Jeans law (note that a = P — 3). 

We see from Figure [T^ that most of the extrapolations 
are not consistent between the spectral bands, with the 
exception of the cut-off PL in 7-rays and the fiattest PL 
in X-rays. Note that a steep PL fit (P — 3) in X-rays 
would need to break sharply above FUV frequencies in 
order to connect the X-ray and FUV spectra. The PL- 
only fit to the Fermi data also appears unlikely, as it 
would need to turn over strongly as it approaches the 
X-ray band. Similarly, the PL fits to the FUV spectrum 
must turn over before the softest energies of the X-ray 
band, but this would be expected for thermal emission. 

It turns out to be possible to fit simple models to the 
whole set of multi-wavelength spectral data. For ex- 
ample, in Figure [131 we have fitted a single cut-off PL 
across the whole range, with a number of thermal com- 
ponents (black-bodies for the UV- X-rays, and the same 
WD atmosphere as before in the IR/optical). Since in 
this case, the WD part is fixed, we only fit for frequen- 
cies > 9 X lO^^^Hz (A < 3300 A). The optical redden- 
ing is fixed E{B — V) = 0.03 and the X-ray extinction 
Nh ~ 7x 10^^ cm~^. We find that the entire pulsar spec- 
trum (excluding WD) is reasonably well fit {x1 ~ 1.2) by 
a single cut-off power-law and three black-bodies. The 
fit parameters are given in Table [6] (see also Figure [T3)) . 
The uncertainty on the coolest temperature is small be- 
cause the that component contributes in both the FUV 
and soft X-rays, whereas in principle additional thermal 
components could be included in the fit; we note that 
the radius for this component is somewhat too small 
for a neutron star. With the gap in the spectrum be- 
tween the FUV and soft X-rays, we are thus not able to 
constrain the surface temperature better than in §3.1.3, 
1.25 X 10^ < T < 3.5 X lO^K. Further data at the soft- 
est X-ray energies (< 0.2 keV) and/or an independent 
measurement of the X-ray extinction would be very use- 
ful. The latter could possibly be obtained from measures 
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Fig. 12. — Multi-wavelength spectral energy distribution of PSR J0437— 4715, showing the raw data (points with error bars). Over-plotter 
are curves showing fits from the different wavebands. For the 7-ray band, we show both a cut-off PL and pure PL fits, with dashed lines 
for the best fits, and the l-cr confidence bounds shown with dotted lines. For the X-rays and optical-UV , we show the range of PL 
components from the various models, where the formal uncertainty in each PL slope is small compared to the differences between models 
(a = —2.5, —3.4, —4 for the optical and F = 1.5, 3 for the X-ray). The thermal (pulsar and WD) components arc not shown for simplicity. 



of ph oto-electric edges in hiKh-resol ution grating spectra 
(e.g., iDurant fc van Kerkwiikll2006h . 

We find that thermal emission from the bulk of the 
neutron star surface, with a temperature of l.SxlO^K 
is supported by the data. Furthermore, we have derived 
a rough temperature distribution on the surface of the 
neutron star, largely cool, but with a much hotter polar 
cap (or caps). In this model, it seems plausible that 
the cut-off PL, which can describe the 7-ray emission 
adequately, is a continuation of the same PL detected in 
the X-rays, and which indeed also makes a contribution 
in the NUV range. 

An alternative to the cut-off PL above could be a 
variable-index PL; e.g., we know from multi-wavelength 
observations of young pulsars th at the spectrum of- 
ten s oftens with increasing energy (jKargaltsev fc Pavlov! 
[2007I) . which can be modeled as 

F^^Nv-^-^^°^''^\ (2) 

where z/23 = (i^/102'^Hz) . This model has the advantage 
of producing a near flat F^, in the optical regime, as is 
seen for non-re cyc led pulsars with higher E but steeper at 
higher energiefj- In our case, we can either require the 
PL component to contribute only in 7-rays and X-rays 

^® An exponential cut-off is still expected in the GeV range, but 
with the steeper model spectrum, the fit is no longer sensitive to 
it for our data. 



by fixing h = 0.05 (Figure [T4|. in which case a = 0.89(2) 
(see Table E]). If we allow h to vary, the best-fit curvature 
is 6 = 0.013(3) and the remainder of the parameters tend 
toward those for the cut-off PL fit. Both of the variable 
index fits give statistically worse than the cut-off PL 
above. 

5. DISCUSSION 

With the amount of MW data gathered and analyzed, 
J0437 has the most comprehensive MW spectrum of an 
MSP to date. Also, we have measured the mid-IR-NUV 
spectrum of the WD companion, where, thanks to the 
precisely determined distance, mass and full spectral en- 
ergy distribution, the WD is one of the best-characterized 
WDs known. 

5.1. WD companion 

The modeling in §3.1 . 21 demonstrates the high accuracy 
of the latest WD atmosphere models. We find that the 
WD has a radius R = (1.87 ± 0.20) x 10'' km and tem- 
perature T = 3950 ± 150 K and a hydrogen surface, for 
a WD of mass M = 0.254 M,^. It is one of the coolest 
WDs observed to date (e.g.. lGates et al.ll2004tlHan et al.] 
i2008: Kihc ct al. 2010|). 

The WD radius is about 40% larger than the theoret- 
ical zero-temper ature, zero-rotation radius expected for 
a pur e He WD (jHamada fc SalpeteiifTOGlt iDriebe et afl 
I1998D . This can be seen as a result of the non-zero tem- 
perature hydrogen envelope and the unknown rotation 
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Fig. 13. — Multi-wavelength fit to the spectral energy distribution of PSR J0437— 4715, with WD atmosphere model, black-bodies (labeled 
BB) and a cut-off power-law (F = 1.563(11), i.'cut = 2.7(5) X 10^^ Hz). The inset shows a zoom of the UV part of the spectrum. 



TABLE 6 

MULTIWAVELENGTH FITS 



Fit 


Tbb (lO^K) 


i?BB (km) 


F 


F^(1023 Hz) 


Other parameters 


xVdof 


CutPL+3BB 
VarPL-f3BB 
VarPL-t-3BB 


3.37(9), 1.70(5),0.34(2) 
3. 36(7), 1.66(4),0. 0285(15) 
3. 35(7), 1.69(4),0. 337(12) 


0.039(4),0.160(13),5.3(8) 
0.040(5),0. 168(10), 76(4) 
0.039(6).0.160(11).5.2(8) 


1.562(13) 
1.893(16) 
1.73(3) 


0.99(18) 
3.9(6) 
3.0(4) 


i/cut = 2.7(6) X 10^3 Hz 
b = 0.05 (fixed) 
b = 0.013(2) 


79/54 
130/55 
93/54 



Note. — Black-body and power-law fit parameters for multi-band models: CutPL (cut-off PL) and VarPL (PL with a variable index, 
see text). Numbers in parentheses indicate l-cr uncertainties in the final digit(s). For the VarPL fits, F = a + 1. 



rate. Very high internal magnetization {B ^ 10^^ G) of 
the WD (where the surface field is only a small fraction 
of the interior field strength) could also result in an ad- 
dition al radius increase of up to 10% (jSuh fc Mathewd 

[2OOOI). 

[Hansen fc PhinnevI (|1998[ ) considered the PSR 
J0437— 4715 system in detail, in order to derive limits 
on the age and initial pulsar spin period (since the 
cessation of accretion). The temperature, surface 
composition and mass of the WD companion to PSR 
J0437— 4715 can be use d to derive an estimate o n its 
age. Using Figure 2 of [Hansen fc PhinnevI (|l998f ). we 
estim ate rwp = 6.0 d = 0.5 G yr. For a 0.242 Mq mass 
WD, iSerenelli et al.l ()200H ) also predicted a surface 
temperature and colors consistent withour findings after 
6 Gyr. 

It is interesting to compare the inferred WD age with 
the pulsar age. A correction must be applied to the ap- 
parent characteristic age of the pulsar, caused by the 



transverse Doppler effect due to the pulsar's transverse 
velocity (i.e., proper motion. IShklovskiilll970( ). For the 
particular case of PSR J0437— 4715, being close to us and 
having a large proper motion, this correction is signifi- 
cant. Specifically, the Shklovskii correction to the char- 
acteristic rate for proper motion, /i, and distance, D, 
is given by {P/P)i = [P / P) — [x^Djc where the sub- 



script i denotes the intrinsic value. Using the precise 
VLBI measurements of \x and D by iDeller et al.l j2008!) , 
we find the corrected spin-down characteristic age is 

- (-^)" = 6.7 ± 0.2 Gyr, assuming 

3 and spin period Pq P — 5.8 ms. 



TSR 



1 



Pi(n-l) 

breaking index n - 

The derived ages are rather similar with one-another, 
surprising given the number of assumptions involved. 
This implies Pq must be substantially faster (at least a 
factor of 1.4 for n = 3) than the current pe riod of 5.8 ms 
(cf., 1.4 ms for the fastest known pulsar; iHessels et al.l 

mm . 
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Fig. 14. — Multi-wavelength fit to the spectral energy distribution of PSR J0437— 4715, with WD atmosphere model, black-bodies (labeled 
BB) and a variable-index PL (a = 0.89, b = 0.05 in Equation 2). The inset shows a zoom of the UV part of the spectrum. 



5.2. Pulsar thermal emission 

We have constrained the temperature of the bulk of the 
pulsar surface to 1.25 x 10^ < T < 3.5 x 10^ K for reason- 
able radii, or 1.5 x 10^ K for the preferred radius of 13 km, 
with about a 15% uncertainty dominated by the uncer- 
tainty in the extinction. It may be possible to further 
constrain this temperature by filling in the spectral gap 
from FUV to soft X-rays, is adequate instrumentation is 
available. 

It is obvious from our MW fits that the thermal com- 
ponent cannot be described well by a single BB or 
NS atmosphere model , sup porting the conclus ions of 
iZavUn fc PaviwI (|T99l ) and iZavlin et al.l (|200l . The 
large width of the thermal peak argues strongly in a favor 
of multi-temperature thermal spectrum, which would be 
the natural consequence of a non-uniformly heated NS 
surface. If modeled as a sum of BBs, the observed spec- 
trum requires a high-temperature (~ 4MK) component 
that appears to be emitted from a relatively small area 
and it also requires a cooler component that must be ra- 
diated from most of the NS surface. Hence, one could 
suggest that the hot component comes from heated po- 
lar caps (see K04 and references therein) while the cool 
component is due to the residual heat stored in the NS 
interior. However, the temperature of this cooler compo- 
nent, ~ 0.15 MK, is surprisingly large for a passively 
cooling NS. This means that either a heating mecha- 
nism operates in the NS interior or the heat deposited 
at the polar caps is efficiently conducted to the rest of 
the NS surface. In modeling the phase-resolved spec- 



trum, iBogdanov et al.l (|2007[ ) suggest some axial asym- 
metry for the hotter components explained by an offset 
dipole magnetic field. 

5.2.1. Heating mechanisms 

After an initial neutrino-dominated cooling epoch, pas- 
sive photon-dominated cooling from the NS surface is ex- 
pected to rapidly reduce the surface temperature to be- 
lo w 10^ K by 1 Myr, an d 10^ K b y 1 Gyr (see the review 
bv lYakovlev fc Pethick|[2003 and lPage et al.l[200l . The 
surface temperature we find is well above this. There 
may be a number of possible mechanisms for internal 
re-heating of the NS which can result in higher surface 
temperatur es of very old puls ars. Such mechanisms are 
discussed in lPage et al.l ()2006[) . but we briefly summarize 
several of the relevant possibilities here. 

• The interior of the NS is thought to contain super- 
fluid, so, as the crust spins down, differential rota- 
tion will be created between the layers, and this can 
release energy by the friction between the 'creeping' 
superfluid vortex lines and the non-superfi uid crust 
to which it is pinned ([Larson fc Linklll99aV bimi- 
lar processes are thought to be responsible for the 
well-known glitching behavior of younger pulsars. 
In late times (t 3> 10® yr), the photon cooling from 
the NS surface balances the frictional heating, and 
the thermal luminosity, independent of the previ- 
ous NS thermal history, is equal to the spin-down 
rate multiplied by the excess angular momentum 
in the core with respect to the crust, L — J\(l\. 
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ILarson fc LinkI (|1999( ) used this to explain the tem- 
peratures of three moderately old (r ^ 10 yr) pul- 
sars. 

• Ohmic heating due to curr ents induced by the NS 
magnetic field (jAguilera et al. 2008). This effect is 
only relevant for high-B NSs. Although the initial 
field may have been much higher than it is today, 
all MSPs have weak surface magnetic field, pre- 
sumably due to processes in the accretion phase. 
Stronger fields may, in principle, be buried in the 
NS interior, so we cannot completely exclude ohmic 
heating. 

• The centrifugal force of rapid rotation affects 
the chemical balance in the NS interior. As 
the NS spin downs, the force decreases, and 
the particle populations are forced into non- 
equilibrium. Chemical reactions restoring this 
equilibrium will relea se heat, a process known as 
rotochemical heating (Reisenegger"1995). For the 
surface temperature for J0437 reported by K04, 
IPetrovich fc Reiseneggeii ()2010f ) developed a spe- 
cific model for heating in this case and concluded 
that rotochemical heating was the preferred heat- 
ing mechanism. Our surface temperature range 
is slightly higher than that given by K04 moving 
J0437 slightly further from the rotochemical model 
cooling curve, but not significan tly enough to alter 
IPetrovich fc Reiseneggeii (|2010D 's conclusions. 



• With the high mass and density of a NS, it 
is possible that dark matter particles could be 
trapped there a nd either decay or annihil ate, re- 
leasing energy (|Bertone fc FairbairnI |2008| ). The 
nature of dark matter particles is currently un- 
known, so the quantitative predictions of such a 
mechanism remain rather uncertain. Dark mat- 
ter, even if it liberates both gravitational and 
rest-mass energy, probably has far too low a den- 
sity in the solar neighborhood to contribute sig- 
nificantly to PSR J0437's 1 0^ K surface emission 
(jKouvaris fc Tinvakovll2010l) . 

Both vortex creep and rotochemical heating mecha- 
ni sms depend on rotationa l kine tic energy. According 
to IGonzalez fc Reiseneggeri (j2010D , only these two mech- 
anisms will be important at time-scales comparable to 
PSR J0437's age. Our measurement of the bulk surface 
temperature, 1.25 x 10^ < T < 3.55 x 10^ K, is more 
precise and slightly higher tha n the temperature found 
by K0 4. In Figure 4, top, of IGonzalez fc Reiseneggeii 
(|201Clf ). this pushed the data point from about 1-cr from 
the upper model curve to about 1.5-ct from the model 
curve. Although the allowed range of surface T is in 
agreement with the more uncertain K04 measurement, 
we stress that the important result of our analysis is a 
much increased confidence in the thermal interpretation 
of the pulsar FUV spectrum (K04 allowed for a non- 
thermal interpretation as well). 

5.2.2. Polar cap heating 

As well as residual heat in the interior and any heat- 
ing mechanism discussed above, the surface temperature 
profile of the pulsar will be affected by the energetics 



of spin-down induced magnetospheric activity. A large 
amount of energy is available from the loss of rotational 
kinetic energy, Ei = 2.8 x lO^'^ergs"^, but the major- 
ity of this will be released as a relativistic wind. Some 
fraction of this power is, however, released as magne- 
tospheric photon emission (which accompany the parti- 
cle cascades). High-energy magnetospheric particles are 
responsible for the non-thermal emission from the pul- 
sar, but it is generally believed that some fraction of 
the particle flux returns to the NS surface along the 
field lines, heating the polar caps. Indeed, this is be- 
lieved to be the cause of soft X-ray thermal emission 
with T > 10^ K a nd proje cted areas much smaller than 
the size of a NS (IHardineil20n"l . In the case of J0437, 
it is reasonable to assume that the hotter black-body 
components fitted above, with a combined luminosity of 
Lhot ~ 10'^" ergs~^, are also polar caps heated in a sim- 
ilar way to other pulsars. The thermal X-ray efficiency 
of the order 77 » 3 x 10""^ is within the ex pected range 
(IGil et al.ll2007tJKargaltsev fc Pavlovll2008l) . 

IHardin i l200W reviews the production of particles and 
radiation in the vicinity of pulsars, and polar cap heating 
from different polar cap emission models. Their mod- 
els for MSPs predict a luminosity similar to that which 
we find for the 'hot' luminosity, Lhot- Given that these 
particles will not penetrate beneath the atomic layers 
(jPatta et al.lll995l ) , we expect the downward particle flux 
to have a negligible effect on the temperature observed 
for the bulk of the NS surface. Heat may, however, propa- 
gate along the surface if the conductivity is high enough, 
and this is suggested by the broadness of the thermal 
spectral peak. According to the flts in Table 7, however, 
the hotpolar cap is small, and the temperatures inter- 
mediate between the hottest areas and the bulk of the 
neutron star surface are emitted from an area still small 
compared to the total. We therefore infer that the con- 
duction along the surface is probably small. 

In short, the surface temperature 1.25 x 10^ < T < 
3.5 X 10^ K is unlikely to be affected by magnetospheric 
processes, and therefore reflects the internal thermal evo- 
lution of the pulsar. Furthermore, the spectrally broad 
thermal emission can be resolved into a function of tem- 
perature versus emitting area (c.f. Figure [T3t . which can 
then be used to constrain the heat conduction in the sur- 
face layers, once such models become available. 

5.3. Magnetospheric emission. 

The multi-wavelength data suggest that a single non- 
thermal component may be responsible for emission from 
the UV to GeV energies. Whether this is indeed a sin- 
gle PL, with a high-energy cut-off, or a variable expo- 
nent would best be determined by filling in the spec- 
tral window between 10 keV and 0.1 GeV, which may be- 
come feasible with upcoming facilities such as NuStar. 
The possible rising PL tail at low energies in Figure [T3] 
raises the intriguing possibility of observing the source 
with ALMA. Unfortunately, MIPS on board Spitzer is 
no longer operational, and the data presented here offer 
no constraint on the PL in the mid-IR. 

The GeV luminosity observed is 8 x 10"^^ergs~^ 
(100 MeV-100 GeV, which contains most of the non- 
thermal luminosity as estimated from Figure [T5|) . This 
gives a magnetospheric radiation efficiency rj — L^/ E ^ 
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0.03. T his is within the ty pical range for Fermi-detected 
pulsars (jAbdo et all 120091 ) . 

The n on-thermal emissio n in GeV appears to be 100% 
pulsed (lAbdo et all IMol ). In the X -ray range, the 
pulsed fraction increases with energy fr-om 30% at 0.3 keV 
to oy er 50% for 2-6 keV (Zavhn 2006; Bog danov et al.l 
I2OOI . This can be explained by a combination of weakly 
pulsed thermal component with a smooth pulse profile, 
and a highly-pulsed magnetospheric component with a 
sharper pulse profile. In detailed modelin g of the e nergy - 
dependent light-curve of J0437 Bogdan ov et al.l (pOOTf ) 
do not explicitly consider the role of the PL compo- 
nent. They do remark, however, that the pulse profile 
for the PL alone must not be too sharp, otherwise its 
effect would also be seen in the light-curves at lower en- 
ergies, where thermal emission dominates. If the PL pho- 
ton index were sufficiently flat (F < 2), then the latter 
restriction could be lifted. 

Under the assumption that the non-thermal emission 
is indeed a single PL component, how plausible is such a 
broadband PL? Assuming that we are seeing synchrotron 
emission from FUV to '^1 GeV, a broad range of particle 
energies and/or magnetic field strengths must be present. 
There are two boundary energies: the energy of PL break 
(~1 GeV) and the lowest energy observed ('^1 eV). With 
primary particles capable of radiating at up to 1 GeV 
energies, a population of lower-energy secondary parti- 
cles will be created, which can emit synchrotron radia- 
tion. The photon energy emitted at the light cylinder 
magnetic field (where the dipolal field is the weakest), 
Bhc = 3 X lO'* G, is approximately Egyn ^ 10~''7^ eV 
(where 7 is the Lorentz factor). Thus a lepton popula- 
tion near the light-cylinder with Lorentz factors in the 
range of a about 100 to a few times 10^ could account for 
power-law emission from optical to GeV, with the slope 
of the electron energy spectrum, p = 2r— lf»2.1. The 
value of the high-energy cut-off is typical for the pulsars 
detec ted by Fermi, either classical or MSPs (e.g.. iBarind 
l2010f). presumably because the B^c is similar for both 
populations. As a more conventional interpretations, the 
7-ray emission could be pro duced by the curvature radi- 
ation or primary el ectrons (|Muslimov fc Hardind [200^ 
iHarding et al.l[2008f l. with the energy of the emitted ra- 
diation dependent on the altitude of the emitting parti- 
cles. H owever, it has been recently pointed out by Baring 
([2010') that under certain conditions, the synchrotron ra- 
diation is expected to dominate. 

If the NUV emission is indeed (at least partly) magne- 
tospheric in origin, then one would expect it to be highly 
pulsed, in contrast to the NS and WD thermal emission 
in the optical and FUV. Detection of pulsations would 
confirm the presence of the non-thermal component over 
many orders of magnitude in energy. Only the Space 
Telescope Imaging Spectrograph (STIS) aboard HST has 
both the required sensitivity and time-resolution for such 
an observation. Detection of pulsation from the ground 
(in the u'-band) is possible, in principle, but very sen- 
sitive instrumentation would be required to disentangle 
the magnetospheric and thermal (both WD and pulsar) 
emission components. Likewise, in the (long- wavelength 
end of) FUV, there is no current instrument but STIS 
capable of time resolution high enough to measure pul- 
sation. 



Whereas non-thermal emission has been observed in 
the optical/NUV spectra of young pulsars such as Vela, 
it has not yet been found in the UV for any other old 
(r > 1 Gyr) pulsar. The late-age reappearance of a mag- 
netospheric component in the optical may in principle re- 
sult from two factors: the lower surface heating rate from 
the downward flux of magnetospheric particles; and that 
the heating of the surface from stored internal heat is ex- 
pected to continue in young neutron stars only to ages of 
~1 Myr. The hot thermal component is therefore much 
weaker in older pulsars, enabling the possible detection, 
as here, of both the cool thermal component from the 
bulk of the surface and NUV non-thermal emission. 

Variable-index power-laws have been used before to 
describe broad-ba nd non-thermal emission in pulsars 
(jGusumanol 120041 ). Young pulsars, such as the Crab 
and Vela, are known to emit non-thermal emission in 
a very broad band. For middle-aged pulsars, extrapo- 
lation of the X-ray PL component s to lower frequencies 
comm only over-predict the optical (|Kargaltsev fc Pavlov! 
120071 ). One might speculate that for those young ener- 
getic pulsars where B^c is higher, the low-energy down- 
turn of the PL falls between the optical and X-rays 
(-E10W oc -Blc)j causing a break in the spectrum. On 
the other hand, poor modeling of the X-ray spectrum 
(where detailed atmospheric models may be more ap- 
propriate than simple blackbody models) may bias the 
non-thermal PL slopes measured. 

When the spectral fluxes of the gamma. X-ray and 
IR-UV data of several pulsars are plotted side-by-side, 
connecting the non-thermal emission with a single power- 
law i n each case does seem plausible (see iDurant et al.l 
I2OIII ). That does not necessarily imply that a single 
component is responsible for 100% of the non-thermal 
emission, but possibly for a large fraction of it. The 
prospect of detecting more pulsars in the optical, and 
improved coverage at hard X-ray energies will certainly 
clarify this picture in the near future. 

6. SUMMARY AND OUTLOOK 

We have obtained broad-band spectral coverage of the 
binary system containing the old MSP J0437— 4715. In 
particular, high-quality optical, UV and X-ray spectra, 
supplemented with broad-band coverage in the IR and 
GeV 7-rays have enabled us to construct the most com- 
plete SED for an MSP to date, as well as accurately 
model the atmosphere of one of the best-observed cool 
WDs. The data enable us to constrain the NS surface 
temperature and hint at a broad magnetospheric com- 
ponent stretching from the NUV to GeV 7-rays. 

Although PSR J0437— 4715 gives an interesting data 
point in the T(r) dependence, the thermal history of pul- 
sars and the dependence of evolution on the mass and 
other parameters can only be found by measuring the 
surface temperatures of other pulsars through UV obser- 
vations. A successful model for pulsar thermal evolution 
would constrain the equation of state of the neutron star 
interior. Similarly, an increasing sample of detailed ob- 
serva tions of WDs with w ell-determined characteristics 
(e.g.. iParsons et ani2011f) will invariably lead to better 
evolutionary models. 

If the extremely broad non-thermal spectral compo- 
nent is real, it ought to be seen by nuSTAR in the 10- 
50keV range. Similarly, ALMA may detect the compo- 
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nent at lowest energies, if there is no low-energy cut-off. 
It will be interesting to see whether similar components 
reveal themselves in the spectra of other millisecond and 
ordinary pulsars when the full multi-wavelength spec- 
trum is analyzed. As more Fermi LAT data accumu- 
lates, more pulsars will be available for multi- wavelength 
follow-up. Multi-wavelength studies appear to be the 
most promising exploration tool in the coming next 
decade. 
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